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Abstract 

We discuss how one can discriminate models with cold dark matter (CDM) and 
baryon isocurvatm'e fluctuations. Although current observations such as cosmic mi- 
crowave background (CMB) can severely constrain the fraction of such isocurvature 
modes in the total density fluctuations, CMB cannot differentiate CDM and baryon 
ones by the shapes of their power spectra. However, the evolution of CDM and 
baryon density fluctuations are different for each model, thus it would be possible 
to discriminate those isocurvature modes by extracting information on the fluctu- 
ations of CDM/baryon itself. We discuss that observations of 21 cm fluctuations 
can in principle differentiate these modes and demonstrate to what extent we can 
distinguish them with future 21 cm surveys. We show that, when the isocurvature 
mode has a large blue-tilted initial spectrum, 21 cm surveys can clearly probe the 
difference. 



1 Introduction 



Current cosmological observations are now very precise to give a stringent constraint on 
the adiabaticity of density fluctuations. However, some contribution of non-adiabatic 
fluctuations, so-called isocurvature fluctuations, is still allowed (see e.g., [1]). Isocurvature 
fluctuations can be generated in cold dark matter (CDM) and baryon sectors, from which 
we can have some insight on the nature of dark matter and/or the mechanism of the 
baryogenesis. This is because isocurvature fluctuations might be produced depending on 
how dark matter/baryon number is generated. Such examples are axion |2HH] and Affleck- 
Dine baryogenesis mechanism |9HT2] in which CDM and baryon isocurvature fluctuations 
can respectively be generatec^J. If such a mode is found to be non-zero, it would give 
invaluable information on CDM/baryogenesis scenarios. 

Currently, stringent constraints on isocurvature fluctuations come from cosmic mi- 
crowave background (CMB) observations such as WMAP [1]. However, CMB cannot 
discriminate CDM and baryon ones since these two isocurvature modes give indistinguish- 
able CMB angular power spectra. The initial non-zero isocurvature mode of CDM/baryon 
affects fluctuations of photons, at which we are looking in CMB observations, only through 
gravity. CDM and baryon can be regarded as the same component as far as the gravity 
is concerned, hence the effects of isocurvature fluctuations arise in exactly the same way 
for CMB anisotropy (except the amplitude which depends on the amount of CDM and 
baryon) . 

Thus even if we see some signature of existence of the isocurvature mode from CMB 
measurements, we cannot tell whether CDM or baryon is responsible for that. From the 
viewpoint of testing dark matter models/baryogenesis mechanism, it is very crucial to dif- 
ferentiate these two modes in some way, which is the issue we are going to discuss in this 
paper. Although CMB measurement cannot do this task, the evolution of density fluctu- 
ations of CDM and baryon themselves are not the same for models with initial CDM and 
baryon isocurvature fluctuations. Thus we may have some possibility for differentiating 
those by looking at (CDM/baryon) fluctuations directly. 

In this paper, we discuss how one can probe the difference between these two isocur- 
vature modes by using observations of fluctuations in redshifted 21 cm line absorption. 
As mentioned above, CDM and baryon fluctuations evolve differently when one assume 
different isocurvature modes (i.e., CDM and baryon ones). Hence, if either one can be 
probed directly in some way, we can in principle see the difference. Since 21 cm fluc- 
tuations effectively probe baryon density fluctuations, the observations of them can be 
utilized for this purposq^. 

Other examples producing CDM/baryon isocurvature fluctuations include the curvaton model [TSHTS] . 
in which the isocurvature fluctuations partially/totally correlated with adiabatic ones depending on how 
CDM/baryon are generated [inillT]. Furthermore, isocurvature fluctuations can be generated in neutrino 
density fluctuations and neutrino velocity |18j . However, we focus on uncorrelated CDM and baryon 
isocurvature modes in this paper. 

The possibility of probing baryon and CDM fluctuations separately using 21 cm fluctuations was 
discussed in [19]. The issue of distinguishing baryon and CDM isocurvature modes with 21 cm surveys 
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The organization of this paper is as follows. In the next section, we discuss the evolu- 
tion and effects of density fluctuations of CDM and baryon in models with isocurvature 
fluctuations. We explain in some detail what aspect differs between CDM and baryon 
isocurvature modes by looking at density fluctuations of CDM and baryon in these modes. 
Then in Section 3, we show that 21 cm fluctuations can in principle differentiate the CDM 
and baryon isocurvature modes, which is impossible with CMB observations. We also 
demonstrate a future constraint on CDM and baryon isocurvature fluctuations by utiliz- 
ing future 21 cm survey and discuss to what extent they can be discriminated. The final 
section is devoted to the summary of this paper. 

2 Evolutions of density fluctuations in CDM and baryon 
isocurvature models 



2.1 Notations 

Before we start the discussion on the evolution of density fluctuations, here we give some 
definitions, in particular, regarding isocurvature fluctuations to set our notation. 

Gauge-invariant isocurvature fluctuations between the components i and j are defined 

as 

5,=-3i7f^-^y (1) 

V pi Pj J 

In the following, we mainly consider isocurvature fluctuations of CDM and baryon with 
respect to radiation, which we write as 

3 3 
5'c = 4 - ^^7, S}, = 5b- -5^. (2) 

Furthermore, in this paper we only consider uncorrelated type of isocurvature fluctuations. 
Thus CDM and baryon isocurvature modes can be respectively identified as the ones with 
non-vanishing Sc = Sc{t — )■ 0) and Sb = Sb{T — )■ 0), with r being the conformal time, at 
the early times when the initial conditions are set for cosmological density fluctuations. 
Then the primordial power spectrum for isocurvature fluctuations is defined as 

VsXm^f5^'\k -k') = —, {S.{k)S,{kr) , (3) 

Ztt^ 

where i = c and b indicating CDM and baryon, respectively, and we parametrize P5. as 

(i) 1 

'PsAk) = VsAko)[^J (4) 



was investigated for the "compensated isocurvature mode" in [20]. See also |21) . 
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with 77, i*-* is the spectral index for the mode i. In this paper, we take the reference scale ko 
as ko = 0.002 Mpc~^. For CDM and baryon isocurvature modes, we denote it as ni^^^'^°'' 
and n^^^^°\ respectively. For the adiabatic (curvature) fluctuations, similarly to the above, 
one usually defines as 

Vdk){2nr6^'\k -k') = ^, Umk'T) , (5) 
where C is the curvature perturbation and the power spectrum Vc_ can also be written as 

Pc(fc)=n(M(^^J , (6) 

with Hs being the spectral index for the adiabatic mode. 

Current cosmological observations such as WMAP indicate that cosmic density fluctu- 
ations are almost adiabatic and pure isocurvature fluctuations of any modes are excluded. 
However, some fraction (contamination) of isocurvature fluctuations are still allowed. To 
characterize the fractions of the isocurvature modes, we define the quantities Va and r;, 



as 



.#3 



' b — 



From current observations, rc is constrained to be less than about 10 % [1]. (The corre- 
sponding constraint on becomes less severe by the factor of (Qc/^b)^-) Thus we use 
Tc = 0.1 as a reference value unless otherwise stated in the following. 



2.2 Evolutions of density fluctuations 



Now we study the evolutions of density fluctuations in models with CDM and baryon 
isocurvature fluctuations and see the difference of those between in these isocurvature 
modes, which could be probed with 21 cm fluctuations. In addition, we also discuss why 
CMB cannot discriminate these two isocurvature modes at linear perturbation level. 

In the following, we investigate density perturbations in the synchronous gauge where 
the metric perturbations are defined as 



ds^ = a\T)[-dT^ + {6ij + hij)dx'dx^], 



with 



hij{x,T) 



dk'^e 



3 Ak x 



kikj - ^6ij ) 6r]{k, r) 



^•^ In some literature such as [Tl[32] , the fraction of isocurvature fluctuations is defined as 

^ VsM 



(9) 
(10) 

(7) 



which is a bit different from the one used here. However, as far as the fraction of the isocurvature 
fluctuations is small, both definitions are almost equivalent. 
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Here r is the conformal time and a is the scale factor. We denote density fluctuations of 
a species i as 5j = Spi/pi. In particular, here we study those of CDM 6c, baryon 6b and 
photons 6.y. We follow the notation of |23j . 

In this gauge, the perturbed Einstein equations are: 

Pr] - = -47rGa^6p, (11) 

= -47rGa\p + P)e, (12) 
h + 2nh - 2k'^7] = -2A7rGa^6P, (13) 

h + 6ri + 2n(h + 6?)) - 2k'^ri = -2A'KGa^{p + P)a, (14) 

where a dot represents the derivative with respect to the conformal time r and 1-1 = a/ a 
is the Hubble parameter defined with the conformal time derivative. 6p and 6P are 
fluctuations of energy density and pressure for the total component. 6' is a variable defined 
as (p + P)6 = J2iiPi + Pi^kvi with Vi being velocity perturbation of species i and a 
represents the anisotropic stress of the total component. The equation for the CDM 
fluctuations is 

5c = -\h. (15) 

Here we note that the velocity perturbation for CDM is set as = so that we can fix 
the gauge freedom in the synchronous gauge. For baryon, we have 

6b = -Ob -hi, db = 'Hdb + Ranc(JT{d^-db), (16) 

where a-r is the cross section for the Thomson scattering and rig is the number density of 
electron, and we also have defined R = Ap^/{3pb)- For photons, the equations are: 

4 2- • /I \ 

6^ = "~2^7 ~ 9j = k^ i -6^ ^ ^7 j ^ aueariO^ — 9b). (17) 

Similarly, for neutrinos, the equations of motion are 

Su = -t^u - h, 9, = e (U, - a^i . (18) 



3 3' \A 

Here we have omitted the equation for the anisotropic stress and higher multipole moments 
for photons and those for neutrinos since they are not so relevant to our discussion and 
we do not discuss in the following. 

During the epoch deep in the radiation dominated era, the time scale of the Thomson 
scattering is much shorter than the Hubble expansion, which dampens anisotropic stress of 
photons and drives the velocity perturbations of baryon and photons the same. Thus in the 
above equations, we can set cx^ = and 9b = 9^ in early times when the initial conditions 
are given. Furthermore, on superhorizon scales, the terms with velocity perturbations can 
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be regarded as smaller by a factor of f i^i . which allows us to neglect those terms at leading 
order of k. Then we obtain the following relations at the leading order in kr: 



h = -26c = -2^6 



(19) 



which shows that the evolutions of 6i follows that of h with some numerical factors except 
a constant term which can be determined by the initial condition for the perturbations. 

On the other hand, inside the horizon, the evolutions of fluctuations are different 
from those on superhorizon scales. Although CDM holds the relation with the metric 
perturbation h with Eq. f[T^ . on sub horizon scales where /cr ^ 1, baryon and photons 
obey the following equations: 



6b + n6b + clk^6b 

H — k'^6^/ 
3 



-h. 

3 



(20) 
(21) 



Since inside the horizon, 6c grows as 6c (x In a during radiation-dominated (RD) era and 
6c (X a during matter-dominated (MD) era, the metric perturbation h also evolves as 
h (X 1/a and a^^^ during RD and MD eras, respectively. Since the scale factor can be 
given by the conformal time r as a oc r and during RD and MD eras, respectively, from 
Eqs. ( l20l) and ( |2T]) . 6b and 5^ show oscillatory behaviors with almost constant amplitudes 
except the epoch of the transition from RD to MD. As will be discussed shortly, these 
behaviors are the same for adiabatic and isocurvature modes except the phases of the 
oscillations due to the difference of the initial conditions of these modes which are given 
when the fluctuations enter the horizon. 

Now we discuss the evolutions of density fluctuations for each mode by showing those 
of CDM, baryon and photons. We start with the adiabatic mode. To identify the adiabatic 
mode, it is convenient to introduce the gravitational (Newtonian) potentials and ip which 
are defined as the metric perturbations in the conformal Newtonian gauge: 



(r) [-dr'^il + 20) + (1 - 2ip)6,jdx'dx^] . 



(22) 



The metric perturbations h and rj in the synchronous gauge can be related to and ip as 

1 



2fc2 



h + Qi] + H{h + Qr]) 



1p = 7] 



(23) 



In the previous section, we have defined the primordial power spectrum for the adiabatic 
mode with the curvature perturbation (, which is related to ip and as 



C = ^ + 



3{l + w) 



n 



(24) 



"^■^ This is not true in the neutrino density and velocity isocurvature modes 
in this paper. 



which we do not consider 
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Figure 1: Evolutions of 5ci Sf, and for adiabatic mode in the synchronous gauge. Here 
we show those for k = 1.0 (top), 0.01 (middle) and 0.0001 Mpc~^ (bottom). 
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^=1 .0 Mpc (CDM isocurvature) ^=1 .0 Mpc (baryon isocurvature) 




Figure 2: Evolutions of (5c, and 5^ for CDM isocurvature mode (left panel) and baryon 
one (right panel) in the synchronous gauge. Evolutions of 5c — Sc and 6c — Sh for CDM and 
baryon isocurvature modes respectively are also plotted. Here we show those for A; = 1.0 
(top), 0.01 (middle) and 0.0001 Mpc"^ (bottom). 
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where w = P/p is the equation of state for the total component. 

The adiabatic mode can be identified with a non-zero value of ( for r — )■ 0, which leads 
to a solution on superhorizon scales as, at leading order in the series of kr |18j . 

VC(0) = (25) 

^/^(O) = 1 5 + 4/ 

"^^ ' 12(15 + 4/^) ' ^ ' 

5c = S, = ^6, = ^6^ = -\k'T\iO), (27) 

9^ = 9, = -1a;V=^C(0), (28) 



36 

4/. + 23 



4/. + 15, 



(29) 



with the fraction of neutrinos in the energy density of radiations being denoted as fi, = 
^ . Here the conformal time r is normalized such that the scale factor is given by a(r) = 
r+r^ with the normalization for a being Oeq = 1/4. Although our main focus in this section 
is the evolutions of density fluctuations in models with isocurvature modes, for reference, 
we show those in the adiabatic mode in Fig. [1] In the figure, cosmological parameters are 
taken as nbh\= Ub) = 0.0227, n^h^{= = 0.108, h = 0.724 and r^eion = 0.089, where 
Qi = Pi(ao)/pcrit is density parameter for a species i, h is the Hubble parameter and Troion 
is the optical depth for reionization. Unless otherwise stated, we assume these values in 
this paper. 

Next we look at the evolutions of density perturbations in the CDM isocurvature mode. 
This mode can be identified with a non-zero value of (5c at r — 0, which we denote as 
Sc- Notice that even in this mode, the evolutions of 6i and h follow Eq. (fT9i) . Then the 
superhorizon solution in early times can be given by [H] 

h/Sc = 4ficr - eny, (30) 

r]/S, = -^fi,r + ^]cr^ (31) 

SJSc = 1 - 2ncT + 3ny, (32) 

Sb = = = i-2n,T + 3ny)Sc, (33) 

9^ = 9t = 9, = -^n.er^S,. (34) 

Notice that the metric perturbations h and t] are induced by the existence of the initial 
isocurvature fluctuations. Importantly, once the metric perturbations are induced, density 
fluctuations are also generated through gravity, which satisfy the adiabatic condition as for 
those of the adiabatic mode. Thus the superhorizon solution of 6c for the CDM isocurvature 
mode has the structure of the form: 

6c = (initial fluctuations) + (induced from metric perturbations). (35) 
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In the left panels of Fig. [21 we show the evolutions of 5c, Sb, as well as 6c — Sc which 
corresponds to the part induced by the metric perturbations in this mode. 

On the other hand, the baryon isocurvature mode can be identified as a mode where 
6b does not vanish in the early times. The superhorizon solution in this mode is 



h/Sb 


= 4Qbr - GQbT^, 


(36) 


V/Sb 


2 2 


(37) 


6b/ Sb 


= 1 - 21]fer + 3^]fer^ 


(38) 




= {-2nbT + 3nbT^)Sb, 


(39) 


6^ = 6b = By 


= -^Qbk^r^Sb, 


(40) 



where the initial value of 6b is denoted as 5";,. The superhorizon solution for baryon 
isocurvature mode above can be obtained by just replacing Qc in those for the CDM one 
with Qb (and Sc with Sb) except that 6b and 6c are replaced with each other. Thus as far 
as 6^ is concerned, its evolutions in CDM and baryon isocurvature modes are the same 
besides some overall factors such as Qc and Qb originating from the amount of CDM and 
baryon (see Fig. [2]) . 

So far we have seen that the initial, superhorizon solutions of all perturbed quantities 
other than 6c and 6b differ only by an overall constant between the CDM and baryon 
isocurvature modes. In fact, this is also true for all subsequent times, which was also 
noted in [21] previously. Here we discuss this issue in some details. 

First of all, the superhorizon (initial) solutions of all perturbation quantities other 
than 6c and 6b are determined by the initial total matter isocurvature perturbation S^ = 
fcSc+fbSb, not by Sc and Sb separately (here fc = ^c/^m, fb = ^b/^m and Qrn = ^c + ^b), 
which can be noticed by rewriting the superhorizon solutions with Sm'- 





h = 






(41) 




V = 








(42) 


6^ 


= 6,= 


g 




(43) 


01 


= e, = - 








(44) 




6m — [l 2r2mT + 3r2m 


T 


21 r- 
J Jm, 


(45) 



where 6m = fc6c + fb6b- (Pure) CDM and baryon isocurvature modes correspond to the 
cases of = and Sc = 0, respectively. Secondly, the time derivatives of the perturbed 
quantities of any fluids are dependent on 6c and 6b only through the metric perturbations 
h and hence the matter perturbation 6m (See Eqs. ( ITSlfTSl) ). Therefore, as far as the 
value of Sm is fixed, regardless of the values of Sc and Sb, the subsequent time evolutions 
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of all perturbed quantities apart from 6b and 6c are the same. For instance, as can be 
seen (except for the overall normalization) in Fig. [21 the evolution of 5^ is the same for 
the CDM and baryon isocurvature modes. This explains that models with the same Sm 
give identical CMB power spectra. Therefore, in the linear cosmological perturbation 
theory, the CDM and baryon isocurvature modes are in principle indistinguishable from 
observations of CMB. Furthermore, they also cannot be distinguished by observations of 
neither the metric nor the total matter perturbations. 

Put in another way, the above can also be understood as follows. From a naive expec- 
tation, one may ask why 6c and 6b would affect the perturbation evolutions in the same 
way even though baryon couples to photons while CDM does not. The reason is twofold. 
First, photons and baryon couple via the elastic Thomson scattering and the momentum 
is just transferred between photons and baryons. Therefore the photon perturbations are 
affected by 9b, but not by 6b (See Eq. ( ITT]) ). Second, both CDM and baryon are pres- 
sureless matter, thus the velocity perturbations of these fluids 9c and 9b are not generated 
from the density perturbations 6c and 6b (See Eq. ( 1TB]) ). which is not true for photons 
or neutrinos (See Eqs. f lTTll and f lTSj) ). These two facts explain that 6b and 6c affect the 
evolution of cosmological perturbations only through gravity in the same manner. 

While the CDM and baryon isocurvature modes cannot be distinguished from the 
perturbations in the photon fluid nor the metric, the evolutions of 6c and 6b are different 
between these two modes as seen from Fig. |2l Although the evolutions of the induced parts 
of 6c and 6b (represented as 6c — Sc and 6b — Sb in Fig. [2]) from the metric perturbations are 
also the same except some constant factor and retain the relation implied by Eq. ( IT9|) . the 
initial nonzero fluctuations of 6c and 6b for CDM and baryon isocurvature modes should be 
superimposed, respectively, on those induced parts. Thus the evolutions of 6c and 6b are 
different, which gives a possibility of differentiating these modes. Although the difference 
between 6c and 6b becomes small on small scales since baryon density fluctuations catch 
up those of dark matter after the recombination and entering the horizon. However, on 
large scales, as seen from the bottom panels of Fig. [21 we can see some differences between 
6c and 6b even at later epoch. Interestingly, we can in principle observe this difference by 
looking at 21 cm fluctuations which effectively probe the fluctuations of baryon. In the 
next section, we discuss this issue in some details. 

3 CDM and baryon isocurvature fluctuations with 21 
cm Survey 

3.1 21 cm fluctuations 

First we briefly summarize the fluctuations in redshifted 21 cm line emission. Even after 
recombination at 2 ~ 1090, the photon and the hydrogen gas are still in kinetic equilibrium 
via the Compton scattering due to a small fraction of free electrons. At around a redshift 
z ~ 300, the interaction becomes ineffective and thereafter the gas cools adiabatically. The 
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gas temperature Tgas scales as a ^, while the photon temperature Ty scales as a ^. At first, 
the 21 cm spin flip driven by the atomic collision is so frequent that the spin temperature 
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tracks Tgas- As gas cools faster than CMB, the atomic collision subsequently becomes 
less effective and the interaction with the abundant CMB photons in turn dominates. 
Then Ts catches up with until the formation of first objects, which is supposed to take 
place at z < 30. Since remains lower than throughout this 'dark age' at redshifts 
30 ^ 2; ^ 300, the redshifted 21 cm line intensity is observed as an absorption compared 
with the CMB intensity at the microwave band [25] . 

The unperturbed brightness temperature of 21 cm line observed now at an observed 
frequency u can be calculated by [251126] 



[l-e 



T - T 
l + z 



where is the redshift of the 21 cm transition at the frequency i.e. I 
with the frequency of 21 cm emission being denoted by vq ~ 1.4 GHz. 
optical depth of 21 cm, which can be given by 



(46) 

1 + z^) = 1^0/ u, 
T2icm{z) is the 



T"21cm(2;) 



z)H 



1.081 X 10"^ xh 



'T 

-7 


-1 








-1/2 


'1 + z 


1/2 




[0.023 J 




L0.I3J 
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(47) 



where Aiq = 2.85 x 10~ s~ is the spontaneous emission coefficient, nui is the number 
density of neutral hydrogen, xhi is the neutral hydrogen fraction, and Um = cUb + ^c- In 
this paper we consider observations at redshifts Zi, > 30 so that complexities arising from 
nonlinear physics can be much reduced. In particular, we consistently set xh = 1 at these 
redshifts, adopting the standard assumption that the onset of the reionization occurs at 
latter epoch. Then Tfe(z/) in Eq. fH6|) can be given by 

1/2 



TbM = 29.50mK 



T —T 



bJb 



0.023J L0.I3J 



-1/2 
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(4^ 



Now we turn our attention to the fiuctuations in Tb- While there can be a number of 
sources generating the 21 cm brightness temperature fiuctuations, the most dominant ones 
are as follows: (1) the monopole source arising from the density fiuctuations of neutral 
hydrogen and relative temperature difference between the spin and CMB, (2) the Doppler 
redshift due to the peculiar velocity of neutral hydrogen gas against the observer. Taking 
into account these two effects, fluctuations in Tb at frequency v observed at position tq 
and conformal time tq in a sight direction h can be given by 



AT. 



T 



■(A, 



A. 



T"21c 



-n 



dv/dr 

n 



Tg determines the number ratio of neutral hydrogens in the triplet and singlet states, ritripiot/'' 
3e~^*/-^% where iJ* is the energy difference between the singlet and triplet states. 



(49) 



dot 
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where Ax = {X — X)/X represents the fractional perturbation of a quantity X in the 
Newtonian gauge. Eq. ( 149|) is formal solution from the line of sight integration, and the 
right hand side should be evaluated at position rv = rQ + hijQ — Ti,) and time t = t^, where 
Tjy is the conformal time at z^. As shown in [27], in a redshift range z > 30, other sources of 
21 cm brightness temperature fluctuations are not very significant (at most percent-level 
at very large angular scales), as long as narrow redshift window functions are adopted. 
Therefore in the following, we adopt Eq. fH9|) as the theoretical prediction of observable 
21 cm fluctuations. 

The observed fluctuations in 21 cm intensity are determined by the fluctuations in the 
spin temperature ATg, as well as A„jjj, At^ and v. Since the rates of the spin flip transition 
due to atom collisions and photon interaction are sufficiently high, the evolution of ATg 
can well be determined by assuming equilibrium. For a comprehensive discussion for the 
evolution equations for ATg and other fluctuations related for 21 cm fluctuations, we refer 
to [27]. 



3.2 3-dim power spectrum from 21 cm tomography observations 

21 cm photons emitted at some redshift z now would be observed at a distinct frequency 
^ = ^o/{^ + z). This allows us to perform tomographic reconstruction of 3-dim power 
spectrum of redshifted 21 cm fluctuations [25 | l26 | [28] . 

Let us consider a tomographic observation of 21 cm fluctuations around a flxed fre- 
quency u and a direction n in the sky. The observed frequency and sky position constitute 
a 3-dim space and we denote its coordinate as (5z/, 0). Once given a background cosmol- 
ogy, there is a one-to-one correspondence between {Su, 0) and the coordinate of 3-dim real 
space {r\\,r±): 

r|| = y{z^)5u, r± = dA{z^)Q, (50) 

where the subscripts || and _L represent the direction parallel and transverse to the line of 
sight, respectively. The conversion factors between (5//, 0) and {r\\,r±) are given by 

In particular, dA{z) is the comoving angular diameter distance. Then a Fourier dual of 
{5u, 0), which we denote as (uy, mj_), can also be related to the wave vector (/cy, A;_l) via 

k\\ = u\\/y{zy), k± = u^/dA{zy). (53) 



From the formal solution Eq. fH9|l . we obtain the power spectrum of 21 cm fluctuations 
observed at around frequency v and sight direction n in A; space 

P^T.ik- V, h) = f 2e-2^— [Poo{k, z,) - 2^i-'Po,{k, z,) + fi^Pn{k, z,)] , (54) 
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where fi is the cosine of angle between n and k, i.e. fi = n ■ k/k. Here Pqq, Pu, Pqi are 
the power spectra of auto-correlations of monopole Aq = ^nui + f'^f i^'^s — ■^^7) and 
dipole moments Ai = -^jf^f^j-f^ source and their cross- correlation, i.e. 

{Ao{k,z)Ao{k',z)) = Poo{k,z){2nf6^'\k + k'), (55) 
{A^{k,z)A,{k\z)) = Pnik,z){27Tf6^^\k + k'), (56) 
{Ao{k,z)A^{k',z)) = Poi{k,z){27rfS^^\k + k'). (57) 

By performing a Fourier transformation of the two-point correlation function of the 21 cm 
fluctuations (ATf,(r)AT;,(r')), the power spectrum defined in the u space, PAniu), can be 
related to one in the wave vector space via 

Note that due to the Doppler effect, the power spectrum of 21 cm fluctuations PAT^iu) is 
not isotropic; it is instead a function of both and u± or, equivalently, both fcy and k±. 

In Fig. [3l we show the 21 cm power spectra for the cases with pure CDM and baryon 
isocurvature fluctuations. To calculate the power spectra for 21 cm fluctuations, we used 
CAMB sources [27j. For reference, we also show expected sensitivity from the Fast Fourier 
Transform Telescope (FFTT) which can probe the redshift z > 30 and is described in [29] . 
An explicit form of noise power spectrum P^^^{u) is given in Eq. fl6Tl) and some detailed 
descriptions will be given in the next subsection. 

Although any pure isocurvature modes are already excluded by current observations of 
CMB, we show pure isocurvature spectra to see to what extent we can tell the difference 
between CDM and baryon isocurvature modes. In the top panels of Fig. [3l we show power 
spectra for the case with = ni^"^°'' = 1. The amplitude of primordial fluctuations 

for CDM mode is fixed as Vs,{ko) = 2.41 x 10"^ at A; = 0.002 Mpc"^ For the baryon 
mode, the amplitude is assumed to give exactly the same CMB angular power spectrum 
as that for CDM (i.e., the initial amplitude is (Qc/^b)^ times large as that for CDM 
mode). The bottom panels show the case of very blue-tilted spectra with n^^^^^^"^ = 
^(biso) _ 2 Although the spectral indices of this kind of value seem somewhat too large, it 
has been claimed that blue-tilted (CDM/baryon) isocurvature spectrum may be favored 
by observations [2211301 - IMj^ . In addition, some axion model producing very blue-tilted 
spectrum has been discussed [12], and in their model, the spectral index for the CDM 
isocurvature perturbations can be as large as 4, depending on parameters. In light of 
these, it would be interesting to investigate isocurvature models with very blue-tilted 
spectra. 

As seen from Fig. [31 the survey of 21 cm fluctuations gives different power spectra for 
CDM and baryon isocurvature fluctuations, which enables us to discriminate these modes. 

A recent analysis puts a constraint on the spectral index as ni'~^^^^^°^ ~ 2.2461q'428 ^M- 
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Figure 3: 21 cm power spectra -Pat6(w||, m±) at z = AO for the cases with pure CDM (sohd 
red hue) and baryon (dahsed green hne) isocurvature modes. The spectral indices are 

1 (CDMiso) (biso) i /; \ j (CDMiso) (biso) o /i NT i 

assumed as ris = nl =1 (top) and ris = nl = o (bottom). In each 

panel, power spectra are shown as function of k± with k\\ being fixed. From left to right, 
M|l is fixed to 0.22MHz~^, 0.72MHz~^ and 2.3MHz~^, which respectively correspond to 
~ 5.4 X 10-3Mpc-\ 1.7 X IQ-^Mpc"^ and 5.6 x lO'^Mpc-^ The size of the baryon 
isocurvature mode is chosen such that the baryon isocurvature mode gives the same CMB 
power spectrum with the CDM one. We also show the FFTT noise power spectra (dotted 
blue line). 
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However, when the primordial spectra are assumed to be scale-invariant (top panels), the 
signals are too weak compared with sensitivities of observations. On the other hand, for 
the case with ni^^^'^°'' = n^^^^°'' = 3 (bottom panels), the amplitudes become large and 
difference in power spectra can larger than the noise power. So we expect that these kinds 
of isocurvature modes can be differentiated with the 21 cm survey of FFTT. 

In fact, even when we consider a realistic case where CDM/baryon isocurvature fluc- 
tuations give subdominant contributions to the total (almost adiabatic) spectrum allowed 
by current observations, we can still probe the difference between CDM and baryon isocur- 
vature modes. In Fig. HI we show the 21 cm spectra for {rc,rb) = (0.1,0) (red line) and 
{rc,rb) = (0,0.1 X (Qc/^b)^) (green line). Contributions from the adiabatic perturbations 
are also included here. The choice of the value of for the latter case is made so that 
these two cases give the indistinguishable CMB angular power spectra. In the figure, we 
also show the expected error in measuring power spectrum 5PATb(w), as well as that for 
the noise P^'^'^iu)- An explicit expression of 5P^Ti,iu) is given in Eq. fl64p . While the noise 
power spectrum is estimated from the detector noise, error accounts for both the cosmic 
variance (assuming a vanilla model with the WMAP 7yr mean parameters [1]) and survey 
setups. The figure indicates that these two cases can be differentiated with the FFTT 
survey. To see this issue more quantitatively, in particular, to what extent we can probe 
the difference with a future 21 cm tomography observation such as FFTT, we make a 
forecast by adopting the Fisher information matrix analysis in the next section. In closing 
this section, we also note that the the difference between the two isocurvature modes are 
more prominent at higher redshifts. This reflects the fact that the density perturbation in 
baryon catches up with one in CDM as redshift decreases, as we mentioned in the previous 
section. 



3.3 Discriminating CDM and baryon isocurvature fluctuations 
with future 21 cm survey 

Now we make a quantitative discussion about the possibility of discriminating CDM and 
baryon isocurvature modes with future observations of 21 cm fluctuations. For this pur- 
pose, we assume FFTT [25] for 21 cm survey. Although 21 cm surveys can probe the 
difference between CDM and baryon isocurvature modes, regarding other cosmological 
parameters, CMB observations are much more powerful. Thus we also include a future 
CMB observation, for which we assume the specification of CMBpol [35] in the Fisher 
matrix analysis. 

Before presenting our results, here we briefly describe the formalism of the Fisher 
matrix analysis. Since we use expected measurements of CMB and 21 cm survey, the 
Fisher matrix Fij can be given by the sum of two contributions: 

^_^(CMB)^^(21cm)^ (59) 

where indices z,j represent cosmological parameters and F^^^^^ and Fj^^^^^^ respectively 
indicate the contributions from CMB and 21 cm observations. We analyze a forecast in a 9 
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k^[Mpc" 



Figure 4: 21 cm power spectra Pat^ , u±) from the mixture of adiabatic and uncorrelated 
matter isocurvature perturbations. From top to bottom, we show power spectra at z = 30, 
40 and 50 in order. In each panel, power spectra are shown as a function of k±, for a 
fixed From left to right. My is fixed to 0.22MHz-\ 0.72MHz-i and 2.3MHz-\ which 
respectively correspond to k\\ ~ 5.4 x 10~^Mpc~^, 1.7 x 10"^Mpc"^ and 5.6 x 10~^Mpc~^. 
Shown are power spectra for the cases with (rc,rfe) = (0.1,0) (dashed green line) and 
{rc,rb) = (0,0.1 X (fic/fife)^) (dotted blue line). The spectral indices for the isocurvature 



fluctuations are assumed as n^J"^^^^"^ = nl"'""' = 3. Notice that these two cases give 
completely the same CMB power spectrum. We also show the FFTT noise power spectra 
(dotted blue line) and expected errors in measuring 21 cm power spectrum (magenta bar). 
Errors are estimated adopting a pixelization, which is logarithmically uniform in {u\\,u±) 
with widths Alnwii ^ AlnM_|_ ^ 1. 
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Frequency (GHz) 


beam FWHM (arcmin) 


AT (/iK arcmin) 


45 


17 


5.85 


70 


11 


2.96 


100 


8 


2.29 


150 


5 


2.21 


220 


3.5 


3.39 



Table 1: Specifications for the mid-cost (EPIC-2m) CMBpoI mission used in the analysis. 
For polarization, AT is reduced by the factor of a/2. We have excluded the highest and 
lowest frequency bands (30 GHz and 340 GHz) for a realistic foreground removal. 



dimensional parameter space {oOc, oob, h, Ug, TVeion, Ag, ni^^^^^°^). The definitions of 
the parameters here are already given in Section I^T^ except Ag, the amplitude of primordial 
fluctuations at the reference scale k = 0.002 Mpc~^. Here we assume that the spectral 
index for baryon isocurvature mode is the same as that for CDM one, i.e., n^^^^"^ = rig^^^"^ 
is adopted in the analysis. 

For the Fisher matrix for CMB Fj^^^^\ we follow the method of [30] and assume 
the projected CMBpol mission [33]. More specifically, we use the specifications for the 
mid-cost (EPIC-2m) mission, which is given in Table [TJ 

The Fisher matrix from a 21 cm tomography survey can be given by [371138]. 

p(2icm) ^ f d^u Ve / dPAni^) \ ( dPAn{u) \ 

J {2nYP'^^^{uY \ d\ )\ dx, )' ^ ' 

where Aj is the i-th cosmological parameter and P^°^ (u) is the power spectrum of the 
total observed 21 cm fluctuations, which is the sum of the signal PaT(,(^) and the noise 
power spectra P^^^f (u), i.e. P'^i^^{u) = PAniu) + P^^f (u). Here Ve = ^fovP is the 
survey volume in G space, where fipov and B are the solid angle in the field of view and 
the frequency band width, respectively. In the analysis, we set fipov = tt and B = 8 MHz, 
regardless of observed frequencies. 

Noise power spectrum of FFTT can be given by 

AT,, l"J - -J?^ f T- 

i FoV J cover ''obs 



= -T7\ 1 f— W"(m±) • (61) 



Here /sky is the observed sky fraction, X^, is the observed radio wave length, Tgys is the 
system temperature, A is the telescope area, and /cover is the fraction of are covered by 
the array. In the analysis, we take /cover = /sky = 1, A = 20km^ and tobs = lyr- W{uj_) is 
the beam response function, which in principle depends on the shape of the telescope. In 
the analysis, we assume a symmetric Gaussian beam and we take 



W{u±) = exp 



X2 



(62) 



17 



Since both Pat,,{u) and P^'^'^{u) do not depend on the azimuthal direction in -u-space, 
we integrate the right hand side in Eq. (!60|) over this direction. In addition, we pixehze 
the u space with widths Au\\ and Au± in the hne of sight and transverse directions, 
respectively. Then Eq. ( I60l) can be approximated as 

p{21cm) _ V- 1 f dPAn{u) \ ( dPAn{u) \ 

-;^/PAT.(nTl dx. )[ dx, J' ^''^ 

where 6Pati,{u) is the error in power spectrum par pixel, which can be given by 



Removal of foregrounds is one of the primary challenges in future observation of 
cosmological 21 cm fluctuations. We assume foregrounds are highly smooth in the fre- 
quency space, and can be removed to sufficient accuracy at > 1/5, where we assume 
Tsys = 200K[(1 + z)/10]'^-^ is dominated by the sky temperature due to the Galactic syn- 
chrotron emission. Note that this choice of lowest hardly affects our results, as we 
focus on a blue-tilted power spectra of isocurvature perturbations. In addition, we also 
exclude modes k > k^a.^ = 0.1 Mpc"^ in the analysis, in order to neglect effects of the 
nonlinear evolution of matter perturbations. 

In Fig. [5l we show projected constraints for and r^, from CMBpol alone (solid red con- 
tours) and CMBpol-|-FFTT (shared regions). Regarding the observed redshifts of FFTT 
survey, we analyzed four different combinations: z = 30 (top left), z = 40 (top right) 
2; = 50 (bottom left) and all these redshifts combined (bottom right). Here the fiducial 
values are assumed as Tc = 0.1 and ri, = and 77,^*-'°'^"^°^ = n^J^^^"^ = 3. Other cosmological 
parameters are assumed as those given in the previous section. Since CMB alone cannot in 
principle make a distinction between CDM and baryon isocurvature models, the constraint 
from CMBpol alone is completely degenerate along the line which gives the same amount 
of isocurvature fiuctuations in term of matter: 



with C = 0.1 X (Q^/Q^)"^. However, as one can see from the figure, by including 21 cm 
surveys, the degeneracy can be removed, which means that we can see the difference be- 
tween CDM and baryon isocurvature modes with 21 cm fluctuation surveys, in particular, 
when isocurvature modes have very blue-tilted spectra such as n^J"^^^^°^ = n^s'^^°^ = 3. 

As discussed in the previous section, the difference between Sc and 6b is more prominent 
at higher redshifts. On the other hand, observational sensitivity is better for lower red- 
shifts. For a flxed band width B and a maximum wave number for the linear perturbation 
evolution /cmax, the increase in the signal overwhelms the increase in noise and the most 
stringent constraint is obtained at the highest redshift z = 50. 
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Figure 5: Expected la and 2(T constraints on and r?, from CMBpol alone (red contours) 
and CMBpol+FFTT (shaded regions). As for FFTT, we show constraints from different 
combinations of observed redshifts in separate panels; 2; = 30 (top left), z = 40 (top right), 
z = 50 (bottom left) and all of these redshifts z = 30, 40 and 50 (bottom right). The 
fiducial values are assumed as (rc,rf,) = (0.1,0) and 77,^*-'°'^^*^°) = 77,^^"^°^ = 3. 
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Figure 6: Same as in Fig. |5]but 
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We also note that the constraint from the combination of three redshifts 2; = 30, 40 
and 50 is at least twice as stringent as one from a single redshift z = 50. This suggest that 
a combination of different redshifts offers to some extent a synergy in discriminating the 
CDM and baryon isocurvature perturbations. This comes from the fact that parameter 
degeneracies between Tc or r^, and other cosmological parameters more or less differ at 
different redshifts. By combining observations at different redshifts, parameter degenera- 
cies are partially removed and this gives severer constraints on Tc and r^. Otherwise, the 
constraint from the combination of three redshifts would be at most as \/3 times tighter 
as one from 2; = 50, which is the most stringent if a single redshift is observed. 

Finally we comment on the dependence of constraints on the spectral indices of pri- 
mordial isocurvature power spectra. In Fig. [6l we show constraints in the r^-rfe plane for 
the case of ni^^'^^^"^ = = 2. As can be expected from the discussion in the previ- 

ous section, constraints on r,, and significantly degrade as the spectral indices becomes 
smaller. For the n^J"^^^^"^ = ni^"^°^ = 2, we may marginally discriminate the CDM and 
baryon isocurvature perturbations for a fiducial model (rc, r;,) = (0.1,0). The significance 
is at around 1 a level from an observation at single redshift 2; = 50. By combining observa- 
tions at all redshifts, the significance can be improved to around 2 a level. However, note 
that to what extent we can discriminate the CDM and baryon isocurvature perturbations 
depends on the maximum wave number for linear evolution. While we adopted a fixed 
^max = 0.1 Mpc~^ at any observed redshifts, the use of larger k^aax may be allowed as red- 
shift becomes higher, which may offer more stringent constraints on Tc and r?,. On the other 
hand, for the case of n^J"^^^^"^ = n^^^^°^ = 1 we found FFTT does not improve constraints 
on Tc and from CMBpol alone regardless of observed redshifts. This is what we expect 
from the top panels in Fig. [3l Thus the CDM and baryon isocurvature perturbations 
would be hardly distinguished if their power spectra are nearly scale invariant. 



4 Summary 

In this paper, we discussed a possibility of differentiating CDM and baryon isocurvature 
fluctuations with 21 cm survey. As is well-known, by using CMB observations, we cannot in 
principle see the difference between those modes although they can give a severe constraint 
on the sum of the size (fraction) of such isocurvature modes. However, once it is confirmed 
that CDM/baryon isocurvature perturbation should contribute to density fluctuations in 
the Universe even if the size is small, it would be very important to differentiate those 
modes since they can give invaluable information on the nature of dark matter/generation 
mechanism of baryon asymmetry of the Universe. 

We showed that 21 cm fluctuation survey, which effectively probes density fluctuations 
of baryon, can in principle differentiate CDM and baryon isocurvature modes as shown in 
Figs. [3] and HI To see this issue in some quantitative manner, we made a Fisher matrix 
analysis and discussed to what extent 21 cm survey can remove the complete degeneracy 
between the fractions of CDM and baryon isocurvature modes which resides in CMB. When 
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isocurvature fluctuations have scale-invariant primordial spectra, it seems very difficult to 
see the difference even with the FFT telescope. However, if isocurvature modes have 
very blue-tilted spectra, which is claimed to be favored by observations p2l[30l - [34] , and is 
predicted in some axion model [12], it would be possible to distinguish CDM and baryon 
isocurvature modes, which we have explicitly shown in Fig. [5l 

The adiabaticity of density fluctuations can give important information on various 
aspects of the Universe. In particular, in light that CDM/baryon isocurvature fluctuations 
might be generated in well-motivated scenarios of CDM and baryogenesis like axion models 
and Affleck-Dine baryogenesis, the differentiation of these two modes, which has been 
shown to be possible by using 21 cm fluctuation survey in this paper, may well give 
important implications for understanding the nature of CDM and baryogenesis scenario. 
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